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Abstract. It is most natural to assume that mysterious Small Ionized and 
Neutral Structures (SINS) in diffuse ISM arise from turbulence. There are two 
obvious problem with such an explanation, however. First of all, it is gener- 
ally believed that at the small scales turbulence should be damped. Second, 
turbulence with Kolmogorov spectrum cannot be the responsible for the SINS. 
\ S Q • We consider, however, effects, that provide spectral index flatter than the Kol- 

fSJ ' mogorov one and allow action at very small scales. These are the shocks that 

arise in high Mach number turbulence and transfer of energy to small scales by 
instabilities in cosmic rays. Our examples indicate that the origin of SINS may 
J> ' be discovered through systematic studies of astrophysical turbulence. 
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Turbulence and SINS 



Various observations covered well in this volume indicate the existence of struc- 
ture in the neutral and ionized ISM on AU spatial scales, causing concern about 
how well we understand the basics of the ISM physics. In the straightforward 
interpretation, these clouds should be extremely overdense and and overpres- 
sured. Evaporation arguments point to very short lifetimes of such clouds, yet 
the AU-scale structures in the diffuse medium seem quite common. In this vol- 
ume the reader can find many articles dealing with the observed properties of 
SINS. Therefore we do not dwell upon this issue. Instead we adopt a broad the- 
^ ■ oretical approach and ask a question whether turbulence can produce enough 

activity on small scales to be considered as a potential source of SINS. 

The attractiveness of turbulence as an origin of SINS is that it can pro- 
duce a generic universal explanation. Indeed, turbulence is really ubiquitous in 
astrophysics, including the ISM and circumstellar regions. However, the most 
popular Kolmogorov turbulence model can not account for SINS. Kolmogorov 
turbulence has a spectrum of E(k) ~ fc~ 5 / 3 , the fluctuations of density at scale 
k ~ 1/7 arisen from advection by turbulence are (d~pk) 2 ~ kE(k) ~ /c~ 2 / 3 , i.e. 
they decrease with the scale. SINS, on the contrary, are connected to rather 
large density contrasts at small scales. This would violate the assumption of 
Kolmogorov turbulence that it is only weakly compressive at large outer scale. 
Fortunately, astrophysical turbulence is not limited to the Kolmogorov spec- 
trum. If the spectrum of density fluctuations scales as E{k) ~ k~ n , n < 1, then 
(dpk) 2 increases with the decrease of turbulence scale. 

As we see from the example above, in terms of turbulence, the origin of 
SINS is a quantitative question related to the spectra of turbulent motions. 
This question cannot be dealt using brute force of computers. Indeed, from 
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the point of view of fluid mechanics astrophysical turbulence is characterized by 
huge Reynolds numbers, Re, which is the inverse ratio of the eddy turnover time 
of a parcel of gas to the time required for viscous forces to slow it appreciably 
For Re ^> 100 we expect gas to be turbulent and this is exactly what we observe 
in HI (for HI Re ~ 10 8 ). In fact, very high astrophysical Re and its magnetic 
counterpart magnetic Reynolds number Rm (that can be as high as Rm ~ 10 16 ) 
present a big problem for numerical simulations that cannot possibly get even 
close to the astrophysically-motivated numbers 1 . The currently available 3D 
simulations can have Re and Rm up to ~ 10 4 . Both scale as the size of the box 
to the first power, while the computational effort increases as the fourth power 
(3 coordinates + time), it is not feasible to resolve the actual interstellar tur- 
bulent cascade to test whether SINS somehow emerge at the end of it. Another 
approach based on (a) establishing scaling relations that can be tested with the 
modern computers and (b) testing these relations with observational data seems 
promissing. In what follows we discuss the part (a). Part (b) is addressed in a 
recent review by Lazarian (2006). 

2. Density spectrum in supersonic MHD turbulence 

A very important insight into the incompressible MHD turbulence by Goldreich 
& Shridhar (1995) (henceforth GS95) has been followed by progress in under- 
standing of compressible MHD turbulence (Lithwick &, Goldreich 2001, Cho & 
Lazarian 2003, henceforth CL03, Cho, Lazarian & Vishniac 2003). In particular, 
simulation in Cho & Lazarian (2003) showed that Alfvenic cascade evolves on 
its own 2 and it exhibits Kolmogorov type scaling (i.e. E ~ /c~ 5 / 3 ) and scale- 

dependent anisotropy of the Goldreich-Shridhar type (i.e. fey ~ k_[ ) even for 
high Mach number turbulence. While slow modes exhibit similar scalings and 
anisotropy, fast modes show isotropy. The density scaling obtained in Cho & 
Lazarian (2003) was somewhat puzzling. At low Mach numbers it was similar 
to slow modes, while it got isotropic for high Mach numbers. 

Density power spectrum shallower than the Kolmogorov one was reported 
in a number of observations (Deshpande, Dwarakanath & Goss 2000, Padoan et 
al 2003) and the relation between the SINs and the shallow power spectrum was 
advocated in Deshpande (2000). 

For high Mach numbers we expect shocks to develop; the density should 
be perturbed most dramatically by those shocks. Naively, one would think 
that density perturbations will have a random-shock spectrum of k~ 2 which 
leaves very little perturbation on small scales. Also, in sub- Alfvenic turbulence 
magnetic field is dynamically important, so we expect to observe some anisotropy 



1 This caused serious concerns that while present codes can produce simulations that resemble 
observations, whether numerical simulations reproduce reality well (see McKee 1999, Shu et 
al. 2004). 

2 The expression proposed and tested in CL03 shows that the coupling of Alfvenic and com- 
pressible modes is appreciable at the injection scale if the injection velocity is comparable with 
the total Mach number of the turbulence, i.e. with {V% + C|) 1//2 , where Va and Cs are the 
Alfven and sound velocities respectively. However, the coupling gets marginal at smaller scales 
as turbulence cascades and turbulent velocities get smaller. 
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Figure 1. Left panel: Mach number is 10, power spectra of: solid - density, 
dashed - velocity, dotted - logarithm of density. Right Panel: Probability den- 
sity function for a density in simulation with Alfvenic Mach number around 
unity and various sonic Mach numbers. 



which is typical for magnetic turbulence. None of these two properties is actually 
observed in simulations. On the contrary, the spectrum was typically rather 
shallow and there was no significant anisotropy. This mystery is resolved as 
follows: shocks in isothermal fluid can have very large density contrasts, up 
to sonic Mach squared, however the conservation of mass and positive sign of 
density allows only small regions with high-density or clumps. These clumps 
totally dominate the spectrum of density. Being close to delta-functions, they 
generate rather shallow spectrum. Being distributed randomly in space, they 
mask any anisotropy originally present due to Alfvenic shearing. 

Our calculations in Beresnyak, Lazarian & Cho (2006, henceforth BLC06) 
demonstrated that the spectrum of density for high Mach number MHD tur- 
bulence is shallow (see Fig. 1). This potentially may be important for SINs, 
although we do not attempt to provide quantitative arguments for this at the 
present stage. The rms deviation of density for a subsonic case is consistent with 
prediction M 2 for low beta (CL03), and the rms deviation of log-density for su- 
personic case is around unity regardless of a Mach number. The distributions 
are notably broader for higher Mach numbers, though. 

Dimension of the high-density structures is between 1 and 2, being viewed 
as a flatted filaments or elongated pancakes. There were no evident preferred 
orientation of these structures along or perpendicular to the mean magnetic field. 
Maximum density value in a Mach 10 data cube was around 3 x 10 3 po- We also 
noted in BLC06 that the randomly distributed high-density clumps suppress 
any anisotropy originating from motions at small scales. 

In magnetically dominated medium that we deal with it is reasonable to 
assume that the corresponding shocks move material along magnetic field lines 
the same way that the slow modes do in subsonic case. The shocks are randomly 
oriented and therefore the clumpy structure that we observe does not reveal any 
noticeable anisotropy. In fact, the density perturbations associated with such 
shocks should not be correlated with the magnetic field strength enhancement, 
which our data analysis confirms. Note, that due to different reasons, this 
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Figure 2. M s <~ 10, Le/i Panel: The isosurfaccs of density, correspond- 
ing to 10 mean densities. Central Panel: Contours of equal correlations 
of log-density, calculated respective to the local magnetic field. They re- 
veal pronounced anisotropy of density structure. Right panel: and r± 
for anisotropic eddies revealed by logarithm of density, and comparison with 
GS95 law, solid M s ~ 10, dashed M s ~ 3. 

correlation also weak for density fluctuations induced by slow modes low Mach 
number MHD turbulence (see CL03). 

Interestingly enough, similar results in terms of shallow density spectrum 
are seen in pure high Mach number hydrodynamic simulations (Kim & Ryu 
2005). This reinforces the idea that the primary mechanism for generating of 
flat spectrum of density is the random multiplication/division of density in su- 
personic flows which, together with mass conservation, form isolated high peaks 
of density. Does this mean that magnetic field is totally unimportant for the 
structure of density in supersonic turbulence? Not at all - we were able to 
reveal the underlying density structure, that comes from Alfvenic shearing by 
using a log-density instead of density for spectra and structure functions. The 
structure function of log-density have shown very familiar Goldreich-Sridhar 
scale-dependent anisotropy, typical for almost any quantity in sub- Alfvenic com- 
pressible turbulence (see Fig. 2). 

Calculations of the density scaling with a more extended sample of simula- 
tions with a variety of orginary and Alfven Mach numbers in Kowal, Lazarian & 
Beresnyak (2006) are consistent with our findings above. The explanation of the 
shallow spectrum of density presented above suggests that if we include gravity 
the clumps will be denser and the power spectrum of density will be flatter or 
could be even rising [5- functions will give the ID power spectrum of k 2 ) 

In terms of the relation to SINs, we may note, that cooling may make 
interstellar gas more pliable to compression than the isothermal gas that we 
used in the simulations. This, could result in more density contrast when the 
original gas is warm. This calls for detailed studies that include cooling. 

3. Alfvenic Slab motions induced by Cosmic Rays 

In this section we describe a new mechanism to transfer turbulent energy to 
small scales bypassing the usual turbulent cascade. This mechanism, described 
in Lazarian & Beresnyak (2006) allow us to have relatively large-amplitude slab- 
type Alfvenic perturbations on small scales. Even though MHD Alfven waves do 
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not perturb density, small-scale kinetic Alfven waves (KAW) do. Also the slab 
Alfven waves in a compressible media is subject to the well-known parametric 
instability that produces density perturbations. 

The close connection between the cosmic ray (CR) power-law and the tur- 
bulence power-law is the point that was stressed by Randy Jokipii in many of 
his presentations (see Jokipii 2001). The physical essence of this relation and, 
in fact, the very understanding of the interaction of CR and turbulence are cur- 
rently the issues of intensive research. It is well known that the propagation 
of CR depends on the scaling of the ISM turbulence. Most of the studies in 
the field use isotropic or slab-like Alfvenic turbulence the origin of which does 
not follow from either theoretical models of MHD turbulence (see Goldreich &; 
Shridhar 1995; review by Cho & Lazarian 2005 and references therein) or numer- 
ical simulations of MHD turbulence (see Cho & Lazarian 2003 and references 
therein). On the contrary, the interactions of Alfvenic turbulence with scale 
dependent anisotropy that is consistent with both theoretical predictions and 
numerics were shown to be inefficient for cosmic ray scattering (Chandran 2000, 
Yan & Lazarian 2002). Using the scalings for fast modes obtained in Cho &i 
Lazarian (2003), Yan & Lazarian (2002) identified fast modes as the principal 
scattering agent for standard models of MHD turbulence, provided that the en- 
ergy is injected at scales larger than a pc. In this model that follows from this 
work the differences in damping of fast modes result in the differences in CR 
propagation in different ISM phases. This induces substantial changes on the 
models of CR propagation in our galaxy (Yan & Lazarian 2004) and in clusters 
of galaxies (Brunetti & Lazarian 2006). 

While the above changes of the picture of CR propagation and possibly 
acceleration are inevitable in the quasi-linear models of scattering by MHD tur- 
bulence injected at large scales, in an astrophysically realistic situation with the 
pressure of CRs close to thermal or magnetic pressure, it is reasonable to ask 
whether the feedback of CRs to MHD turbulence is important. This question 
was addressed in Lazarian Sz Beresnyak (2006). CRs react very differently from 
ordinary gas when the magnetized fluid is compressed, if the scale of compres- 
sion is less than CR mean free path. Such scale compressions or expansions 
change only the component of CR momentum that is perpendicular to magnetic 
field (due to the so-called "adiabatic invariant" conservation), while the compo- 
nent of the momentum parallel to magnetic field stays the same. Such state is, 
however, unstable. 

We calculated the rate of CR kinetic instability that transfers energy to 
small scales, creating small-scale Alfvenic perturbations that are not a part of 
the global MHD cascade. These perturbations are more like waves moving along 
magnetic field lines and thus are very different from highly anisotropic GS95 
Alfvenic modes. 

Quantitatively, we consider a power-law distribution of CRs Fq ~ p~ a ~ 2 
where a is conveniently defined as the power-law index for a one-dimensional 
distribution (or particle density). For example, around the Earth a ~ 2.6 up to 
the energies of 10 14 eV. The the growth rate of the cosmic-ray- Alfven gyroreso- 
nance instability (henceforth GI) can be estimated as: 

7CR(fc||) = ±w pi —AQ, (1) 
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where ucr(p > mQ/ku) is the number density of CRs with momentum larger 
than the minimal resonant momentum for a wave vector value of fey, m is the 
proton mass, n is the density of the thermal plasma, uj p i is the ion plasma 
frequency. Q is a numerical factor that depends on the index of cosmic rays 
a. The ± sign corresponds to the two MHD modes. We shall concentrate on 
the Alfven mode, corresponding to the plus sign, as those are less subjected 
to linear damping. We show that when anisotropy is created by compressive 
turbulence, the anisotropy factor A = (p± —p\\)/p\\ is small and changes its sign 
on the scale of the mean free path, depending on two competitive mechanisms - 
scattering which tends to isotropize momentum distribution, and magnetic field 
compression which tends to make it oblate or prolate. 

Assuming that anisotropy arises from the turbulent compressions with am- 
plitude Sv at large scale, the factor A is equal to 2Sv/va, where va is the Alfven 
velocity. Then the expression for the instability rate can be written as 



8v 



where r p is a Larmor radius of a CR resonant with a particular wave vector 
fey = mCl/p, ro is the 1 GeV proton Larmor radius and 

7 1 / B \ f 4- 1(T 10 cm- 3 \ 

This CR instability gets energy from external turbulent compressions of CR at 
the mean free path scale and directly transfers it to scale of the CR Larmor 
radius. 

There are two non-linear processes limit the growth rate of the instability. 
First of all, the magnetic perturbations generated at the Larmor radius of CR 
backreact on CR by limiting the mean free path A. If the change of magnetic field 
direction is <\> ~ SB/B the scattering that is a random walk requires N ~ l/cf 2 
interaction and 

A~ Nr p ~r p /<f? ~r p B 2 /(SB) 2 , (4) 

where we designated SB as the magnetic field perturbation pertaining to a par- 
ticular wavenumber, i.e. SB 2 « E(k)k. We can consider SB as a function of 
either k or the resonant Larmor radius r p (see Longair 1994). SB grows as the 
instability grows, which in turn reduces A. On the other hand, it is the mean 
free path A which determines the scale at which compressions of the magnetic 
field are important. This can be understood as follows: the CR distribution 
"remembers" the perturbed value of the magnetic field and its anisotropy only 
during the time the typical particle travels its mean free path. Once particles 
scatter significantly, the anisotropy of the distribution is effectively "reset" . As 
a result only low amplitude motions on scales less than A excite the instability 
which limits the degree of anisotropy A attainable. 

The instability grows as d(8B 2 )/dt = Jcr^B 2 ) where the injection of en- 
ergy is happening at the scale of the mean free path, i.e. 
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Figure 3. Left Panel: Decorrelation of a plane, k± — Alfven wave by 
turbulence. Lower picture shows the energy density of a wave in cylindrical 
k-space. Alfven waves were injected at = 17. Wave energy is being trans- 
ferred in the direction of k± axis, which is typical for decorrelation by MHD 
turbulence. Upper plot shows decay rate of the wave vs its wavenumber. Right 
Panel: Energy density of compressive modes and Alfvenic slab-type waves, 
induced by CRs. The energy is transferred from the mean free path scale 
to the CR Larmor radius scale. If the mean free path falls below compres- 
sive motions cutoff or feedback suppression scale, the spectrum of slab waves 
becomes steeper. See more on feedback in Lazarian & Beresnyak (2006). 



where eqs. (J2J) and (@J) were used. We see that according to the above equations 
5B perturbations will grow as i 1 / 2 ^ thus reducing A virtually to r p . If not for 
other processes, it is possible to show that this suppression would reduce A to 
virtually the CR gyroradius, inducing Bohm-type diffusion of CRs. 

However, the collection of waves with different wavelengths created by the 
instability is also subject to steepening with the rate of 7steep 

-(5B/B) 2 k\\VA, 

which combined with the effect of limiting A above provides: 



8B rj /2 (T p \ (M-a+2)/(2 M+ 2) 



B ~ l V(2m+2) l m/(2^+2) y rQ J 



(6) 



For a = 2.6 and \i = 1/3 Eq. © produces a rather shallow spectrum of pertur- 
bations, E(k) sa {5B) 2 /k ~ k~ - 8 . 

The calculations above are valid provided that A is larger than the com- 
pressive mode cutoff scale / cut . If, on the other hand, l cut > A the compression 
for the instability is supplied from the eddies at the damping scale, namely, 
bvjvA ~ ('cut/£) 1 ^ 3 (A/7cut)- In this case instead of eq. © one gets 

5B ( rV 2 \ V4 fr \ ( 3 " a )/ 4 



b V^ 1/4 i' l/4 ^ 1 7 At)/4 / 

which, for the same value of a = 2.6, corresponds to a steeper spectrum of 
E(k) ~ £r L2 . 

The ambient Alfvenic turbulence provides yet another source of non-linear 
damping of the waves generated by the instability. This process is analogous 
to the suppression of the streaming instability (Yan & Lazarian 2002, Farmer 
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& Goldreich 2004). Our numerical testings of the scalings of the corresponding 
damping rate, i.e. 

7t urb ~ -k ± v ± ~ -VAkTL- 1 '* ~ -VAr-WL- 1 ' 2 , (8) 

are shown in Figure 3 (left). The encouraging agreement of Eq. |H] with the 
results of MHD simulations makes us confident in using this rate, which up 
to the sign convention, coincides with the prediction in Farmer & Goldreich 
(2004). Combining eqs. © and (jHJ) we find that for a > 5/3 our instability 
is damped for all scales, larger than 

r^^ro^-^Lrf^-^. (9) 

Therefore the spectrum of plane Alfven waves given by eq. ((BJ will protrude 
from r PjCr i t down to r Pim i n which corresponds to minimum energies of CRs. 



4. Summary 

Astrophysicists always expect to see Kolmogorov turbulence. These expecta- 
tions may not be met in reality. We presented two examples of astrophysical 
processes that provide distinctly non-Kolmogorov spectra. One case is the den- 
sity in supersonic MHD turbulence, another is the velocity and magnetic field 
perturbations at small scales. This work is suggestive that a possible explanation 
of the SINS should be sought invoking non-Kolmogorov turbulence. 
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